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ABSTRACT 

We measure the intrinsic relation between velocity dispersion (cr) and luminos¬ 
ity {L) for massive, luminous red galaxies (LRGs) at redshift z ^ 0.55. We achieve 
unprecedented precision by using a sample of 600,000 galaxies with spectra from the 
Baryon Oscillation Spectroscopic Survey (BOSS) of the third Sloan Digital Sky Survey 
(SDSS-III), covering a range of stellar masses M, > We deconvolve the effects 

of photometric errors, limited spectroscopic signal-to-noise ratio, and red-blue galaxy 
confusion using a novel hierarchical Bayesian formalism that is generally applicable to 
any combination of photometric and spectroscopic observables. For an L-cr relation of 
the form L oc cr^, we find /3 = 7.8 ±1.1 for cr corrected to the effective radius, and a 
very small intrinsic scatter of s = 0.047 ± 0.004 in logjQ a at fixed L. No significant 
redshift evolution is found for these parameters. The evolution of the zero-point within 
the redshift range considered is consistent with the passive evolution of a galaxy pop¬ 
ulation that formed at redshift z = 2 — 3, assuming single stellar populations. An 
analysis of previously reported results seems to indicate that the passively-evolved 
high-mass L-cr relation at z ^ 0.55 is consistent w ith the one measured at z = 0.1. 
Our results, in combination with those presented in iMontero-Dorta et al.l (|2014D . pro¬ 
vide a detailed description of the high-mass end of the red sequence (RS) at z ~ 0.55. 
This characterization, in the light of previous literature, suggest that the high-mass 
RS distribution corresponds to the “core” elliptical population. 

Key words: surveys - galaxies: evolution - galaxies: kinematics and dynamics - 
galaxies: statistics - methods: analytical - methods: statistical 


1 INTRODUCTION 

In the 60’s and 70’s, several empirical scaling relations be¬ 
tween the kinematic and photometric properties of early- 
type galaxie s (ETGs) were ident i fied. Than ks to the sem¬ 
inal w ork of IPiorgovski fc DavisI lll987l) and IPressler et all 
digsiti . today we know that these relations are different pro¬ 
jections of the so-called fundamental plane (FP) of ETGs, 
which is a thin plane outlined by the occupation of ETGs in 
the three-dimensional space spanned by velocity dispersion, 

* E-mail: amontero@astro.utah.edu 


effective radius and surface brightness, i.e., logjQ cr, logj^Q Re, 
and logjQ(7}e, respectively. 

Of particular interest among these scaling relations is 
the L-cr relation, a specific two-dimensional projection of the 
FP that relates the luminosity L and the central stellar ve- 
loc ity dispersion cr o f ETGs. This relation was first reported 
by Minkow^ d 19621). us ing a sample of only 13 ETGs, and 
by Morton fc ChevalieJ dl973l ) a decade later, although no 
quantification was provided. The relation was first quant i- 
fied in a sample of 25 galaxies bv iFaber fc Jacksorj dl976l) . 
as a power law in the form L oc cr^, and has been commonly 
called the Faber-Jackson Relation (F-J relation) since then. 
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As a link between a distance-independent quantity a and an 
intrinsic property L, the F-J relation immediately became 
a useful distan ce estimator and conse q uently a cosmologi¬ 
cal p robe ('e.g.. Ide Vaucouleurslll982al la: [^aturel fc GarnieJ 
Il99:j) . 

With the emergence of large-scale structure (LSS) 
galaxy surveys in the last decades, the size of the ETG 
samples available increased dramatically, especially at low 
redshift, i.e. ^ ^ 0.1 (mostly with the Sloan Digital Sky 
Server, SDSS. lYork et ani2000l ). It was then confirmed with 
statistical significance, that the L-cr relation at intermediate 
masses/luminosities approximately follows a canonical F-J 
relation with a slope of ~ 4 (note that the words “slope” an 
“exponent” are commonly used interchangeably, as the re¬ 
lation is often expressed in the form M oc logj^g a, where M 
is the absolute magnit ude). This result is reported in, e.g., 
ISernardi et al.l (l2003al ) and IPesroches et al.| ll2007l) . along 
with a typical scatter in logjg cr of ~ 0.1 dex. 

Almost since the very discovery of the L-a rela¬ 
tion, however, it became clear that the slope depended 
on the luminosity range of the sample under analy¬ 
sis. A value of > 4.5 has been measured in luminous 
ETGs (e.g.. ISch echterl jl980l: iMalumuth fc Kirshn^ Il98ll : 


ICappellari et ah 20131 : ^^onnendv^Bende^ 20W) and 
2 in faint ET Gs fe.g.. iTo nrvI Il98ll : iDavies et al.l 19831 : 
Held et ^Il992l: Ide Riicke et al.ll2005l : iMatkovic fc GuzmanI 


20oJ7"A~confirmation of the curvature of the L-ct relation 

towards the high-mass end with high statis tical significance 
has been possible thanks to th e SPSS (see |Desroches_et^l] 
200?!^ Ih Yrje_fc Bemardl l2009al : iNigoche-Netro et al.l 201^ 


Bemardi et al. I I 2 OIII) . The z ~ 0.1 sample of iBernardi et al.l 
20 111. as an example, contains ~ 18,000 ETGs with 


logiQ M* > 11.2/11.5. As the F-J relation, where L oc ct"*, is 
therefore a particular, canonical case that only appears to 
hold for morphologically selected samples at intermediate 
mass ranges, in the remainder of this paper we adopt the 
generic “L-a relation” terminology. 

The curvature of the L-cr relation is a consequence of 
the curvature of the FP itself. In fact, the phenomenology 
of the FP has turned out to be rather complex. Enough evi¬ 
dence has been gathered that its characteristics depend, not 
only on luminosity/stellar mass but, to a greater or lesser 
degree, on a variety of other galaxy properties. In addition, 
the sensitivity of the measurements to selection effects and 
low-number statistics has often lead to contradictory re¬ 
sults. In terms of redshift evolution, it appears clear that 
the zero-point of the FP must evolve from z ~ 1 in a way 
that approxi mates that of a passively-evolving galaxy pop - 
ulation (e.g., iBernardi et al.lf2003bl : Ijorgensen et al.l 20 od) . 
However, there are some indications tha t the FP might 
be si g nihcantly steeper at higher redshifts llJorgensen et al.l 
I 2 OO 6 I : iFritz etld] l2009l ). whereas for the intrinsic scat- 
ter, conclusions are so fa r somewhat contradictory (e.g., 
iNigoche-Netro et al.l 1201 ll measure a decrease in the scat¬ 
ter o f the the L-cr relation, whereas results from IShu et ahl 
I 2 OI 2 I appear to indicate the opposite). It seems that the 


properties of the FP also depend, to some extent, on the 
wavelength range probed, which is an indication that stellar 
population properties have an impact on the FP. While these 
pro perties appear to rema i n fairly unchanged in th e opti¬ 
cal llBernarcli et al.l l2003bl : iHvde fc BernardU l2009lJ) , there 
are clear indicati ons that the FP is intrinsically different in 
the infrared (e.g., Pahre et al.lll99^ : iLa Barbera et aI]l201Cll : 


iMagoulas et al.N2012h . The dependence of the FP on envi¬ 

ronment has also been the subject of a number of works 
(see e.g. Ti;:eii_et_al, 2001 Bernaxdi_et^ 2003b[;^ Redn^t ^ 


|2004| . 1 2005 1 _ Demcol6_et_aJj 2005 jLa Barbera et al.ll2010l . for 

the FP: lFocardi fc Malavasill2012l for the L-cr relation). This 
is not by any means a complete list; a number of other galaxy 
properties have been investigated in recent years. 


While the use of the FP and its projections as a cosmo¬ 
logical probe has been largely eclipsed by other techniques, 
the FP has gained increasing attention in the last decades 
as a source of observational constraints for galaxy formation 
and evolution theories. In this sense, the variation and evolu¬ 
tion trends with respect to galaxy properties in the slope and 
intrinsic scatter of the FP and it projections are believed to 
be the imprints of non-homological physical processes that 
occur during the formation and evolution of galaxies. Under¬ 
standing these observed trends provides crucial insight into 
these fundamental processes. In this sense, recent cosmo¬ 
logical simulations have investigated the impacts of various 
physical processes such as major/minor mergers and disc in¬ 
stabilities on shaping the i nternal structure and kinematic 
prop e rties of ETGs (e .g.. lOser et al.| 2 0121 : IShankar et al.l 
I 2 OI 3 I : IPosti et al.ll2014 : IPorter et al.1 2014 ). 


In this paper, we use data fr om the Baryon Osci lla- 
tion Spectrosc opic Survey (BOSS. IPawson et al. I I 2 OI 3 I) of 
the SDSS-HI (lEisenstein et ahl 1201 ll ) to measure, for the 
hrst time, the high-mass end of the L-ct relation at 0.5 < 
z < 0.7. This is a continu ation of the work presented in 
iMontero-Dorta et al.1 ll2014h . hereafter MD2014, where the 
intrinsic colour-colour-magnitude red sequence (RS) distri¬ 
bution is deconvolved from photometric errors and selection 
effects in order to compute the evolution of the RS luminos¬ 
ity function (LF). An important conclusion of MD2014 is 
that, at hxed apparent magnitude, and for a narrow redshift 
slice, the RS is an extremely narrow distribution (< 0.05 
mag), consistent with a single point in the optical colour- 
colour plane. This work is intended to measure the L-ct rela¬ 
tion that this photometrically distinct population obeys. At 
intermediate redshifts, the BOSS capability to characterize 
the massive RS population is unrivaled, with a huge sample 
of more than 1 million luminous red galaxies (LRGs) with 
stellar masses M* > 10 ^^Mq. No other previous survey or 
sample has been able to probe this population with compa¬ 
rable statistics, hence the unique value of the measurements 
reported here (find a preliminary analysis from BOSS in 
IShu et aLll2012l . where incompleteness is, however, only par¬ 
tially addressed). Importantly, our intrinsic RS was identi- 
hed in MD2014 using exclusively photometric information, 
i.e. our red-blue deconvolution is based on the phenomenol- 
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ogy of the colour-colour plane, and not on any morphological 
classihcation. We will, therefore, use the “RS” terminology 
instead of the “ETG” terminology when referring to our 
sample/results. 

The mass range covered by BOSS has been hard to 
probe at z ~ 0.55. In fact, most of the information that 
we have about the high-mass end of the L-cr relation comes 
from the SDSS at 2 ~ 0.1 or from small samples of very- 
nearby ETGs, using high-resolution observations (the latter 
being strongly affected by low-number statistics and selec¬ 
tion effects). One of the most important discoveries from 
these studies is that the slope of the L-cr relation is steeper 
at higher masses/luminosities, as mentioned above. In recent 
years, a picture that attempts to explain this mass depen¬ 
dence has emerged. The curvature in the scaling relations 
has been associated with a characteristic stellar mass scale 
of ~ 2 X 10^^ Mq. This scale, which was first reported at 
high significance by iBernardi et al.l (1201 il l with the SDSS, 
is thought to be relat ed to a change in the assembly history 
of galaxies (recently, IBernardi et al.l |2014| have shown that 
this scale is also special for the late-type galaxy population). 


High-resolution imaging has shown that the high-mass 
scale marks the separation between two distinct ETG 
populations: core ellipticals are defined by the fact that 
the central light profile is a shallow power law separated 
by a break from the outer, steep Sersic function profile, 
whereas in coreless ellipticals (also known as “power-law” 
or “cusp” elliptica ls) this feature is not present (see e.g. 
iLauer et al.l l2007bl lal. It has been shown in small sam¬ 
ples that core ellipticals dominate at the high-mass end. 


(iKormendv & Bender! 

19961: iFaber et alj 1 19971: iHvde et alJ 

I 2 OO 8 I: CaoDellari et al 

20131: Kormendv & Bender 2013f). 


Importantly, this bimodality in the central surface bright¬ 
ness profile extends to a variety of other properties. 
The distinct characteristics of each of these types, in¬ 
cluding the fact that core ellipticals obey an L-a rela¬ 
tion with a significantly st eeper slope (iLauer et al ] |2007bl: 
iKormendv fc Benderl l2013l i. have been associated with 2 
different evolutionary paths for these objects. Core ellipti- 


less mergers (iDesroches et al. 2007 

: von der Linden et al. 

2 OO 7 I: iHvde et alJ l2008l: ILauer et al. 

2007al: IBernardi et al. 

2 OIII: CaoDellari et al.l l2013l: Kormendv & Benderl 201311. 


whereas coreless ellipticals might have undergone more re¬ 
cent episodes of star formation llKormendv et all l2009l re¬ 
view evidence that they are formed in wet mergers with 
starbursts). 


To complement the work done in MD2014, here we 
present a novel method to combine photometric and 
spectroscopic quantities in low signal-to-noise (SN) large- 
photometric-error samples that we call Photometric De- 
convolution of Spectroscopic Observables (hereafter, PDSO). 
The PDSO is a hierarchical Bayesian statistical method that 
allows us to combine the velocity disp ersion likelihood func¬ 
tion measurements of IShu et al.l ll2012l i with the photometric 


red/blue deconvolution of MD2014 to provide the most pre¬ 
cise measurement ever performed of the high-mass end of the 
intrinsic L-cr relation within the redshift range 0.5 < 2 < 0.7. 

This paper is organized as follows. Section provides an 
overview of methods and motivations. In Section|3]we briefly 
describe the target selection for the galaxy sample that we 
use, the BOSS CMASS sample (13.111 . and the computa- 
tion of stellar v elocity dispersion likelihood functions from 
IShu et ciD (l2012l i (13.211 . Section |4] is devoted to summarizing 
the results of MD2014 regarding the intrinsic RS colour- 
colour-magnitude distribution. In Section [5] we present our 
PDSO method n a general form {531 and we addresses 
the application of our method to the BOSS CMASS sam¬ 
ple (15.211 . In Section [6l we describe our aperture correction 
procedure. In Section [3 we present the best-£t parameters 
for the a - apparent magnitude relation (IZ3I, discuss the 
effect of addressing completeness and the red/blue popula¬ 
tion deconvolution (IL^ and present our the L-a relation 
results (IZ3J. In Section |S] we compare our measurements 
with previous results from the literature (El) and discuss 
on the physical implications of our measurements (18.211 . Ei- 
nally, we summarize our main conclusions and discuss future 
applications in Section[9] Throughout this paper we adopt a 
cosmology with Dm = 0.274, Da = 0.726 and Do = lOO/i km 
s“^ Mpc“^ with h = 0.7( 1 fWMAP7. iK omats u et al.ll201ll L 
and use AB magnitudes (lOke fc Gunn 1983h . 


2 OVERVIEW OF METHODS AND 

MOTIVATIONS 

The statistical power of BOSS to cover the very-massive RS 
population is unrivaled at 2 ~ 0.55, with a sa mple of ^ 1 
millio n galaxies (in the latest data release, see lAlam et al.l 
I 2 OI 5 II with stellar masses > 10 ^^Mq an d a median stel¬ 
lar m ass of M* ~ (as m easured bv iMaraston et ^ 

I 2 OI 3 I . assuming a lKroup^ { 2 OOII I initial mass function). The 
samples, however, present significant challenges, including 
low SN ratio for the spectra, large photometric errors and a 
selection scheme that allows for a fraction of bluer objects 
that increases with redshift. The photometric issues are ad¬ 
dressed in MD2014, where we photometrically deconvolve 
the intrinsic red sequence distribution from photometric er¬ 
rors and selections effects. Our red-blue population deconvo¬ 
lution allows us to characterize completeness in the CMASS 
sample, which is the main LRG sample from BOSS, cover¬ 
ing a redshift range 0.4 < 2 < 0.7 (see next section for a 
complete description of the sample). This characterization 
allows us to analyze the luminosity function and colour evo¬ 
lution of the LRG population. 

The aforementioned analysis in MD2014 is performed 
within the framework of a Bayesian hierarchal statistical 
method that is aimed at constraining distributions of galaxy 
properties, instead of individual-galaxy properties (since an 
object-by-object approach is discouraged by the character¬ 
istics of the BOSS data). The same philosophy is applied in 
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this work to compute the high-mass end of the L-cr relation. 
The main steps of our analysis are: 

• Development of a general Bayesian hierarchal statistical 
method that combines the photometric red-blue deconvolu¬ 
tion and selection function from MD2014 with probability 
density information from a spectroscopic observable. We call 
this formalism photometric deconvolution of spectroscopic 
observables or PDSO. Our method is aimed at constraining 
the hyper-parameters of a model for the joint pdf of sur¬ 
vey galaxies in physical parameter space, by marginalizing 
over the physical parameter likelihood functions of individ¬ 
ual galaxies given the survey data. 

• Application of the above formalism to the computation 
of the best-£t intrinsic L-a relation from BOSS. The idea is 
to parametrize the intrinsic distribution in L-cr space, and 
use the PDSO method to constrain the parameters of this 
distribution: the slope, zero-point and the intrinsic scatter 
of the L-cr relation. 

• Evaluation of the redshift evolution of the best-£t pa¬ 
rameters that define the L-cr relation within a suitable red- 
shift range, namely 0.5 < 2 < 0.7. 


The results of the above analysis will add important 
constraints to the evolution of massive RS galaxies. In 
addition, the PDSO method will lay the foundations for 
future BOSS studies, where the intrinsic distributions of 
spectroscopically-derived quantities can be determined. In 
the broader picture, an important goal of this paper is 
to complement the detailed characterization of the main 
statistical properties of the LRG population initiated in 
MD2014 that will be eventually used, in combination with 
N-body numerical simulations, to investigate the intrin¬ 
sic clustering properties of these systems and the halo- 
galaxy connection in a fully consistent way. The connec¬ 
tion between galaxies and halos will be performed by apply¬ 
ing the techniques of halo occupation distributio ns (HOD: 
e.g.. lBerlind fc Weinberg|[2002l : IZehayi e t al.ll200 5l~) arid halo 
abundance matching (HA M: e.g.. Vale &: Qstrikeil |2004 
iTruiillo-Gomez et al.l [ 2 OI 1 I L These future applications are 
addressed in more detail in the last section of the paper. 


3 THE DATA 

3.1 The CMASS sample 

In this work we make use of both spectroscopic and pho¬ 
tometric data^^om the Tenth Data Release of the SDSS 
(DRIO, lAhn et ^ l2014fl . which corresponds to the third 
data release of the SDSS-III program and the second re¬ 
lease that includes BOSS data. We choose the DRIO in- 
stead of the rec ently published Data Release 12 (DR12, 
lAlam et al.l 1201^1 . in order to be consistent with the lumi¬ 
nosity function results shown in MD2014. The spectroscopic 
DRIO BOSS sample contains a total of 927, 844 galaxy spec¬ 
tra and 535, 995 quasar spectra (this is a growt h of almost 
a factor two as compared to the SDSS DR9, lAhn et ^ 


l2012h . The baseline imaging sample is the final SDSS imag¬ 
ing data set, which contains, not only the new SDSS-III 
imaging, but also the previous SDSS-I and H imaging data. 
This imagining da ta set was released as part of the DR8 
dAihara et al.|[20llh . These imagining programs provide five- 
band ugriz imaging over 7600 sq deg in the Northern Galac¬ 
tic Hemisphere and ~ 3100 sq deg in the Southern Galac¬ 
tic Hemisphere. The typical magnitude corresponding to 
the 50% completeness limit for detection of point sources 
is r = 22.5. The following papers provide comprehensive 
information about te chnical aspects of the SDSS survey: 
iFukugita et al.l ll l996ll desc ribes the SDSS u q riz ph otometric 
system; Gunn et al. 1 1998l l an d I Gunn et al.[ ll2006h describe 
the SDSS c amera and the SDSS telescope, respectively; 
ISmee et al.l (l2013ll provides detailed information about the 
SDSS/BOSS spectrographs. 


The catalog that we used to compute the RS LF in 
MD2014 is the DRIO Large Scale Structure catalog (DRIO 
LSS). The DRIO LSS, which is basically built from the 
BOSS spectroscopic ca talog and is thoroughly described in 
lAnderson et al.l 120141 1. contains a small number of galax¬ 
ies from the SDSS Legacy Survey. The SDSS Legacy Survey 
includes the SDSS-I survey and a small fraction of the SDSS- 
II survey. In the previous SDSS programs, the spectra were 
observed through 3 arcsec fibers, while the aperture size is 
only 2 arcsec for BOSS. This varying aperture size makes the 
LSS catalog slightly heterogeneous as far as the velocity- 
dispersion measurement is concerned. To avoid this prob¬ 
lem, we opt to use the spectroscopic catalog, which is very 
similar to the LSS catalog, and essentially maps the same 
intrinsic population as that described by the RS LF pre¬ 
sented in MD2014. We have also checked that SDSS Legacy 
galaxies are confined to a very small region of low redshift 
and high luminosity within the redshift-luminosity space. 
Importantly, they are mostly found at z < 0.5, which is 
below the redshift range where our LF results are reliable 
(0.52 <z< 0.65, see MD2014). 


We restrict our analysis to the GMASS (for “Gonstant 
MASS”) spectroscopic sample of the BOSS spectroscopic 
catalog. This sample is built within the official SDSS-III 
pipeline by first applying an imaging pre-selection to en¬ 
sure that only detections passing a particular set of qual¬ 
ity criteria are chosen as targets. Secondly, a set of colour- 
magnitude cuts is applied to the resulting catalog, intended 
to select the LRG sample required to effectively measure 
the BAO within a nominal redshift range 0.43 < z < 0.70 
(the sample extends slightly beyond these limits). In a sim¬ 
ilar way, the low-redshift LOWZ sample (not used in this 
paper) is constructed, covering a nominal redshift range 
0.15 < 2 < 0 .43. For more informat io n on the BOSS selec - 
tion refer to lEisenst ein et ah] (l201lll . iDawson et al.l (l2013fl 
and iReid et al.l I 2015[ i (also a summary is provided in 
MD2014). Importantly, the GMASS selection allows for a 
fraction of ~ 37% of blue cloud objects, as measured by 
MD2014 . The stellar masses fo r the red population, as mea¬ 
sured bv iMaraston et ^ (l2013fl . are M, > peaking 
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at M* ~ as suminR a Kroupa initial stellar mass 

function llKroupall200ll 'l. All colours quoted in this paper are 
model colours and all magnitudes cmodel magnitudes. 

The total number of unique CMASS galaxies with a 
good redshift estimate and with model and cmodel apparent 
magnitudes and photometric errors in all g,r and i bands 
in the catalog is 549,005. The mean redshift in the sample 
is 0.532 and the standard deviation 0.128; approximately 
~ 7.5% and ~ 4.5% of galaxies lie below and above the 
nominal low redshift and high redshift limits, i.e., z = 0.43 
and z = 0.70, respectively. For a complete discussion on the 
selection effects affecting the CMASS data, see MD2014. 

3.2 Velocity dispersion likelihood functions 

One of the key ingredients in our L-cr relation study is the 
likelihood function of the central stellar velocity dispersion. 
The method for d etermining this l ikelihood function is de¬ 
scribed in detail in IShu et al.l ll2012tl . Here we provide a brief 
summary. For every galaxy in our sample, the line-of-sight 
stellar velocity dispersion within the central circular region 
of radius 1 arcsec, is measured spectroscopically by fitting a 
linear combination of broadened stellar eigenspectra to the 
observed galaxy spectrum (note that the typical seeing for 
BOSS is 1.5 arcsec). Thus, for the galaxy, the of the fit 
as a function of the trial velocity dispersion is converted into 
the likelihood function of velocity dispersion with respect to 
the observational data di as 

p(di|logiocr) oc exp[-x?(logioCr)/2]. (1) 

The best-fit velocity dispersion and its uncertainty can then 
be inferred from the x^ function. 

Figure [T] illustrates the type of velocity dispersion data 
that we have in BOSS. In the bottom panel, the best-fit cen¬ 
tral velocity dispersion and its uncertainty are shown in a 
scatter plot for all red CMASS galaxies within the redshift 
slice z = 0.55 ±0.005. Here, a simple colour cut g — i> 2.35 
is used to appro ximately remove blu e galaxies in observed 
space ffollowing [Masters et al.ll201ll l. This is just illustra¬ 
tive, as in MD2014 we demonstrate that a simple colour cut 
is not efficient in terms of isolating the intrinsic RS, as large 
photometric errors scatter objects in and out of the colour 
demarcation. At z = 0.55 the colour cut removes 20% of 
the sample. The fraction of intrinsically non-RS objects in 
the CMASS sample at the same redshift according to the 
population deconvolution of MD2014 is 38%. 

The two quantities shown in Figure [T] are inferred from 
the velocity-dispersion likelihood functions. Figure [T] shows 
that the distribution is centered around logj^g cr ~ 2.35, or 
220 km/s, and A logj^g cr ~ 0.06 dex, or 30 km/s. More than 
two-thirds of the subsample have velocity dispersions be¬ 
tween 120 km/s and 300 km/s with 20 — 50 km/s uncer¬ 
tainties. The top panel displays an example of a typical x^ 
function of a galaxy with both velocity dispersion and un¬ 
certainty consistent with the aforementioned central values. 

As the broad shape of the 'xf function displayed in the 
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Figure 1. An illustration of the BOSS velocity dispersion data. 
Top panel: a typical curve as a function of the trial cr 

(D.O.F.=1736). This particular galaxy within the redshift range 
2 : = 0.55 di 0.005 is chosen to have best-fit velocity dispersion 
and uncertainty matching the densest region shown in the bot¬ 
tom panel plot. Bottom panel: scatter plot of the best-fit velocity 
dispersion in logarithmic scale and its uncertainty (also in loga¬ 
rithmic scale, black dots) for red CMASS galaxies {g — i > 2.35) 
within the redshift range 2 : = 0.55 di 0.005. Contours that enclose 
30%, 67%, 95%, and 99.7% of this subsample are overplotted in 
red. 

top panel of Fi gure ^ indicates (and as also emphasized in 
IShu et aI]l2012^ ■ a point estimate for the velocity dispersion 
is only partially informative, even after a dedicated treat¬ 
ment is adopted for the new yf calculation. Note that an 
appreciable fraction of this subsample have best-fit velocity 
dispersions equal to 0 km/s due to the limitations of fitting 
low SN spectra. Using the likelihood function propagates all 
information to the higher level analysis of the entire popu¬ 
lation. That is one of the main motivations for employing a 
hierarchical Bayesian statistical approach in this paper. 


4 BIMODALITY IN THE COLOUR-COLOUR 
PLANE: THE INTRINSIC RED SEQUENCE 
DISTRIBUTION 

A key element in the computation of the intrinsic L-cr rela¬ 
tion for the RS from the CMASS sample is the underlying 
intrinsic RS magnitude and colour distribution. This aspect 
is thoroughly addressed in MD2014, where we present an an¬ 
alytical method for deconvolving the observed (g-r) colour 
- (r-i) colour - (i)-band magnitude CMASS distributions 



log a [km/s] 

z=0.55 ± 0.005 
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from the blurring effect produced by photometric errors and 
selection cuts. The CMASS sample comprises a consider¬ 
able fraction of bluer objects, that can scatter into the red 
side of the colour-colour plane due to photometric errors. In 
MD2014, this aspect is treated by modeling the BC and the 
RS separately (red/blue deconvolution), which allows us to 
correct the RS intrinsic distribution from the contamination 
caused by BC objects. Importantly, this modeling is not per¬ 
formed from the basis of previous assumptions about “blue” 
or “red” objects based on stellar population synthesis mod¬ 
els, but is intended to describe the bimodality found in the 
colour-colour plane. With this consideration, these compo¬ 
nents present the following characteristics: 

• Red Sequence (RS): The RS is so narrow that is consis¬ 
tent, within the errors, at fixed magnitude and for a narrow 
redshift slice, with a delta function in the colour-colour plane 
(width < 0.05 mag), with only a shallow colour-magnitude 
relation shifting the location of this point. The results re¬ 
ported in this paper regarding the L-a relation are based on 
this intrinsic distribution. 

• Blue Cloud (BC): The BC is defined as a background 
distribution that contains everything not belonging to the 
RS and is well described by a more extended 2-D Gaussian 
in the colour-colour plane. The RS is superimposed upon 
the BC, that extends through the red side of the colour- 
colour plane. Again, the name “BC” is not meant to imply 
that this distribution contains only blue, young objects; the 
BC is a spectroscopically and photometrically heterogeneous 
population to which other types of ETCs can pertain, such 
as dusty ellipticals not belonging to the narrow RS. 

The intrinsic distribntion of magnitudes for the RS is 
the key ingredient in the computation of the RS L-a rela¬ 
tion, and this is provided in MD2014 in the form of the RS 
LF. For the sake of convenience, and given the very nar¬ 
row redshift slices that we consider, we will constrain the 
logj^g a — mi (apparent magnitude in the i band) relation 
(note that this is not the L-a relation, which involves abso¬ 
lute magnitudes). As we show in the following section, for 
constraining this relation, it suffices to know the shape of 
the intrinsic distribution of apparent magnitudes, which is 
given by a Schechter Function of the form: 

nsch(mi,( 7 i»,m*,a}) = 0.41og(10)(/>* x 

( 2 ) 

where a = —1 and 0* = unity, for the sake of simplicity. 

The assumption that a = — 1 is dictated by the narrow 
magnitude range, that prohibits fitting for a. The intrinsic 
RS i-band magnitude distribution, as a function of redshift, 
can be obtained by inserting the following linear relation for 
the redshift evolution of the characteristic magnitude m* 
into Equation [2] 

mf®(z) = (4.425 ± 0.125) (z - 0.55) + (20.370 ± 0.007) (3) 


Similarly, for the BC: 

= (4.011 ± 0.178) (z - 0.55) + (20.730 ± 0.010) (4) 

All the relations provided in this section come from the 
analysis presented in MD2014, although they are not explic¬ 
itly reported there. By using these linear relations instead 
of the best fit values at each redshift we avoid introducing 
unnecessary noise into the analysis. These relations are ob¬ 
tained within the redshift range 0.525 < z < 0.65, where 
selection effects are less severe. Note that the BC LF was 
not reported in MD2014 due to the extreme incompleteness 
affecting the BC in the CMASS sample. Instead, it is used as 
a means of accounting for the contribution of the BC in the 
sample and of subtracting this contribution from the RS LF. 
Following the same strategy, we compute the logj^g a — mi 
relation for the BC (and later the BC L-cr relation) only as 
a means of subtracting the contribution of the BC on the 
RS logjQ a — mi relation. 

The last ingredient in the computation of the RS 
logj^g a — mi relation, as we show in the next section, is 
the fraction of BC objects in the CMASS sample, i.e., fuue, 
which is also provided in MD2014. Again, to avoid intro¬ 
ducing noise in the computation, we use the following linear 
relation for the redshift dependence of fbiue'- 

fbiue{z) = (0.890 ± 0.051) [z - 0.55) -b (0.381 ± 0.003) (5) 


5 PDSO: FORMALISM FOR THE 

PHOTOMETRIC DECONVOLUTION OF 
SPECTROSCOPIC OBSERVABLES 

5.1 The method 

The L-cr relation is a relation between a spectroscopic ob¬ 
servable, the stellar velocity dispersion a, and a photomet¬ 
ric observable, the luminosity L (or, in a more practical 
fashion, the logj q cr and the absolute magnitude M). As 
previously mentioned, measuring velocity dispersions from 
BOSS galaxies is hindered by low SN spectra, which dictates 
that instead of point measurements we use likelihood func¬ 
tions for the velocity dispersions. The photometric aspect 
of the L-cr relation computation presents challenges as well: 
as discussed in MD2014, the observed CMASS distribution 
is strongly affected by photometric errors and selection ef¬ 
fects. In this section we present a formalism that combines 
our likelihood measurements for the velocity dispersions and 
our red/blue deconvolution of photometric quantities. This 
formalism, which we call PDSO, can be used for the pho¬ 
tometric deconvolution of any spectroscopic observables, so 
we will present it in a general way. 

Let us start by assuming that we have a sample of galax¬ 
ies indexed by i, each of which has a spectroscopic data vec¬ 
tor di (in our case, velocity dispersions) and a photometric 
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data vector Ci (colours and/or magnitudes). These galax¬ 
ies are selected (i.e., for measuring the spectra) according 
to some photometric cuts described by P{c), which is the 
probability (between 0 and 1) of selecting a galaxy with 
photometric data vector c. P(c) in the case of the BOSS 
CMASS sample can be straig htforwardly derived from the 
CMASS selection scheme (see iDawson et al.llioi^ b 

Assume that we have quantified the probabilistic map¬ 
ping from intrinsic (noise-free) photometric values C into 
observed (noisy) photometric values c, which we denote as 
p(c|C). In MD2014 this mapping is obtained by modeling 
the covariance matrix of magnitudes and colours using SDSS 
Stripe 82 multi-epoch data. 

Assume that these galaxies have been quantified in 
terms of a number of one or more photometrically distinct 
population components indexed by k, with the results ex¬ 
pressed as the intrinsic number density of galaxies per unit 
C in component k, denoted by UkiC) such that nk{C) dC 
has units of spatial number density. In reality, as we will 
show below, this quantity can also have units of number. 
The above characterization can come in the form of the de¬ 
composition into RS and BC galaxy colour-luminosity func¬ 
tions. Note that nk{C) need not be normalized over C, and 
indeed may be divergent. 

As a result of the previous analysis, we have derived 
the individual probability density functions (PDF’s) of each 
component in observed space within the selected sample, as 
well as the fraction of sampled objects in each component. 
To express these in terms of previously introduced quanti¬ 
ties, we introduce the following function for notational con¬ 
venience: 

<l>kic) = P{c) J dCp{c\C)nk{C). (6) 

This is the number (density) of galaxies in observed photo¬ 
metric space for component k. 

The number of galaxies in the sample in component k 
is 

Nk= j dc^k{c) (7) 

The fraction of sample galaxies pertaining to compo¬ 
nent k is 

fk=Nk Y. ■ ( 8 ) 

/ k! 

Note that fk would correspond to fuue (and 1 — fbiue) 
within our deconvolution framework. The PDF of galaxies 
in observed photometric space in component k is 

Pfe(c) = <E>fc(c)/7Vfc. (9) 

The full PDF of the observed sample in observed pho¬ 

tometric space is 

P(c) = ^/fcPi:(c). (10) 

k 

A derived function that we will make use of below is the 


PDF of a sample galaxy in intrinsic photometric space given 
its position in observed photometric space, which we can 
derive by noting that the joint PDF of observed and intrinsic 
photometric vectors for component k is proportional as 

Pk{C,c) o:P(c)p{c\C)nUC) (11) 


so that then 


Pfe(C|c) 


Pk{C,c) 

JdCpk{C,c) 

P(c)p(c|C) nfc(C) 

/ dCP(c)p(c|C)nfe(C) 
p(c|C)nfc(C) 

J dCp{c\C)nk{C) 


( 12 ) 

(13) 

(14) 


Importantly, pk{C\c), which is the function that we will 
use ultimately, is independent of some of the choices that we 
make about the intrinsic distributions, nk{C). In particular, 
any normalization will cancel off, so we can simply use the 
intrinsic number counts shown in Equation [2] where volume 
is not taken into account and />* = 1 unity by dehnition. 

Now assume that there is some parameter or vector of 
parameters x that can be measured from the spectroscopic 
data vector d; in the sense that we can write down and 
compute the function 


p(di|xi) 


(15) 


for each galaxy i. In our case, x is the velocity dispersion, a 
and p(di|xi) is proportional to the function exp[—x^(cr)]. 

Next we assume a parameterized model for the varia¬ 
tion of the spectroscopic observable within the sample popu¬ 
lations as a function of photometric values. This is expressed 
as 


p*,(x|C;t). (16) 

The vector t will denote the “hyperparameters” that de¬ 
scribe this PDF. Our goal is to infer the elements of t. To 
do this, we proceed to express the likelihood function of t 
given the spectroscopic data and the photometric data. In 
our framework, the hyperparameters t only affect the prob¬ 
abilities of the spectroscopic observables given the photo¬ 
metric observations. 


r(t| {di},{Ci}) : 

= p({di}|{ci},t) 

(17) 

]^p(di|ci,t) 

i 


(18) 

/ dxp(di > 

i 

:)p(x|ci,t) 

(19) 

/ dxp(di > 

i 

k 

(20) 

/ dxp(di > 

i 

/ rfCpfc(x|C,t)pfe(C|c7 
k 

1 (21) 


At this point, we have arrived at an expression on the 
right-hand side entirely in terms of quantities that we have 
introduced above, and we can proceed to map and/or max¬ 
imize the likelihood function of the hyperparameters t. 
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5.2 Application to the computation of the L-a 
relation in BOSS 

For the sake of convenience, we proceed by applying our 
formalism to the computation of the logj^Q a — rtii relation 
from the CMASS sample. Our spectroscopic observable is 
therefore the logarithm of the velocity dispersion, i.e. logj^Q cr, 
and the intrinsic distributions (n*,) are represented by the 
Schechter number counts of Equation with a redshift de¬ 
pendence given by Equations O and 2] for the RS and the 
BC, respectively. This choice implies that C and c corre¬ 
spond to the i-band apparent magnitude, rm, in intrinsic 
and observed space, respectively. 

The parameterized model for the variation of the 
spectroscopic observable (logj^g cr) within the sample pop¬ 
ulations as a function of photometric values (rrii) en- 
codes the login (j — mj relation. Motivated by results from 
iBernardi et al.l ll2003ah . we approximate the intrinsic dis¬ 
tribution of velocity dispersions at fixed L as a Gaussian 
distribution in logjg cr with mean < logjg cr > and intrinsic 
scatter s. For component k this has the form: 


-(logio cr- < logipcr >fc)^ 

( 22 ) 

The mean of the velocity dispersion for component k, 
i.e.,< logj^g cr >*,, is assumed to follow a linear relation with 
apparent magnitude, rm, of the form: 


Pfc(logiocr|mi;tk) = exp 

V 


blurring produced by an average seeing of 1.8 arcsec has 
been assumed. 

In order to perform a realistic AC we need to take into 
account the variation of as a function of apparent magni¬ 
tude, rrii, for each redshift slice, which may affect the slope 
of the logj^Q a — rrii relation. To this end, we fit a linear re¬ 
lation to the mean observed i-band logj^g Re measured by 
the pipeline (in arcsec) as a function of rrii El- For a given 
redshift slice, this relation takes the form: 

< logiQ Re >= a + b{mi — 19) (25) 

where we have shifted the reference magnitude to rrii = 19, 
similarly to Equation 1231 for consistency. Figure [2] displays 
the observed i-band log^g Re in arcsec as a function of rrii for 
the redshift slice centered at z = 0.55, in contours enclosing 
67%, 95% and 99.7% of the entire sample, respectively. The 
squares show the mean values in magnitude bins of 0.1 mag 
and the errors the l-cr scatter around the mean. The solid 
line shows the linear fit to the mean values that we use to 
correct, on average, our logjg a — rrii relation. 

By using Equation 1241 it can be easily demonstrated 
that the aperture-corrected ci parameter, i.e. cf^, is related 
to parameter a in Equation [25] in the following way: 


0 ?^= = Cl-0.048a-logio (0.98) (26) 

And, similarly, for cj'^: 


< logio ^ >k= ci,k -f 2.5 -I- C2,k{mi - 19) (23) 

This expression, as we will show in following sections, 
can be easily transformed, within our framework, into the 
L-(j relation, which is expressed in terms of absolute magni¬ 
tudes. 


6 APERTURE CORRECTION 

BOSS velocity dispersions are measured within the 2 arcsec 
diameter aperture of the BOSS fibers. As we move to higher 
redshift within the CMASS sample, the angular size of the 
hber probes progressively larger physical scales. This effect 
is accounted for a posteriori, by applying an aperture correc¬ 
tion (AC) to the best-ht relations obtained by maximizing 
the likelihood function of Equation 1211 By assuming a de 
Vaucouleurs proHle for the variation of the surface bright¬ 
ness as a function of apparent distance to the center of the 
galaxy, we have obtained the following relation: 

(robs/a{< Re) = 0.98 (i?e/Raperture)° °'*® (24) 

that relates the observed velocity dispersion aobs that we 
measure in BOSS, the velocity dispersion averaged within 
the effective radius. Re, and the effective radius itself, in 
arcsec. As part of the derivation of the above relation, the 


= C2 - 0.0485 (27) 

Within the same redshift slice, more luminous galaxies 
are larger in size (see Figure [2]), which implies 5 < 0. The 
AC correction by dehnition tends, therefore, to steepen the 
L-cr relation. 

In order to avoid introducing any extra noise we use 
a linear fit to the values of a and b as a function redshift, 
a{z), b{z). This is shown in Figure 2] for the redshift range 
of interest, 0.5 < z < 0.7. The linear fit that we obtain for 
a{z) is: 

a{z) = (1.47 ± 0.16)z -b (-0.61 ± 0.09) (28) 

and for b{z): 

b{z) = (-0.57 ± 0.32)2 -b (-0.08 ± 0.19) (29) 

Note that the error on the slope of b{z) is large, so we 
will consider also the case where b{z) is constant and equal 
to the mean value within the redshift range 0.52 < z < 0.65, 
i.e. b(z) = -0.402 ± 0.020. 

^ iBeifiori et al.|[2014l present an analysis on the evolution of the 
effective radius in BOSS, but the Re - magnitude relation neces¬ 
sary to derive the aperture correction is not reported. 
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Figure 2. The observed i-band logjQ Re in arcsec as a function of 
the i-band apparent magnitude mi for the redshift slice centered 
at z = 0.55, in contours enclosing 67%, 95% and 99.7% of the 
entire subsample, respectively. The squares show the mean values 
in magnitude bins of 0.1 mag and the errors the l-cr scatter around 
the mean. The solid line shows a linear fit to the mean values. 


A typical value of b{z) — —0.402 implies a correction 
on the slope of the logj^Q a-rm relation (the C 2 parameter) 
of ~ 0.02. Within our framework, the slope of the logjQ cr- 
rrii relation at a given redshift slice coincides with the slope 
of the L-cr relation. As a reference, typical values for this 
slope are ~ —0.1 (or equivalently 4 for the exponent {3 in 
the form L oc cr^; this is the F-J relation). This implies 
that the AC has a significant effect on the slope of the L- 
cr relation when computed from BOSS. The main idea is 
that brighter (and hence larger) galaxies have their velocity 
dispersions corrected by a different factor than fainter (and 
hence smaller) galaxies, for a given fixed angular aperture, 
which affects the slope of the L-cr relation. 

It is interesting to compare the effect of adopting a dif¬ 
ferent aperture correction on the logj^Q a — rm relation. Equa¬ 
tions M and [271 in combination with the values of parame¬ 
ters a{z) and b{z), dictate the sensitivity of the AC to the 
exponent in Equation 1241 By varying this exponent between 
a reasonable range, we can evaluate the effect on the zero- 
point and slope of the log^g a — rrii relation. We have chosen 
a range between 0.04 (corresponding to the value derived by 
Jor gensen et aLlll995ll and 0.066 (found bv ICappellari et al.l 
2006h . The majority o f values adopted in t he literature are 
within this range (e.g. iMehlert et al .1120031 estimate a value 
of 0.06). 

Using Equations 1261 and given that the typical varia¬ 
tion of a(z) within the redshift range considered is ~ 0.2 
(see Figure (3)), we find that the variation in the AC for the 
zero-point, AAC, within the redshift range considered (note 
that a is a function of redshift), for an exponent of 0.048, 
is 0.0096 dex, in absolute values. Increasing the exponent of 



Figure 3. The redshift dependence of parameters b (upper panel, 
squares) and a (lower panel, squares), representing, respectively, 
the slope and the zero-point of the mean relation between the ob¬ 
served i-band logj^g Re (in arcsec) and the i-band apparent mag¬ 
nitude mi. Errors correspond to uncertainty in the fit associated 
to each parameter. The solid lines show the linear fits that we use 
to correct our L-cr relation measurements. 


the AC function to 0.066 would translate into an increase 
in AAC, for the zero-point, of a 38%. On the other hand, 
decreasing the exponent of the AC function to 0.04 would 
result in a AAC 17% smaller. Although the net effect is 
small, these variations can modify the zero-point - redshift 
trend. In Section Es] we discuss the redshift evolution of the 
zero-point, concluding that the effect on the zero-point of 
the uncertainties on the AC is significant, given the narrow 
redshift range that we probe and the mild evolution that we 
measure for the zero-point. 

With regard to the slope, from Equation [57] and Fig¬ 
ure [3] as mentioned above, we find that the typical AC on 
this parameter is ~ 0.02 (i.e. b{z) x 0.048), basically in¬ 
dependent of redshift. Adopting a range of values for the 
exponent of the AC correction 0.04 — 0.066 would translate 
into corrections within the range 0.016 — 0.026. Although 
the effect is not negligible, adopting a different AC would 
not modify the main conclusion of the paper, in terms of 
the steep slope of the L-cr relation, in any qualitative way 
(see following sections). 


7 RESULTS 

7.1 Best-fit parameters for the logjg a-rrii relation 

By maximizing the likelihood function of Equation 1211 we 
obtain the best-fit values for the hyperparameters t. These 
parameters are the following, for each component k (RS and 
BC): Cl,*;, the zero-point of the logj^g a-rrii relation, that cor¬ 
responds to the mean logj^g cr at mi = 19; C 2 ,fc, the slope of 
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this linear relation, and s^, the intrinsic scatter. The opti¬ 
mization of the likelihood function has been performed from 
2 = 0.40 to 2 = 0.70, using a bin size of A 2 = 0.01. This 
redshift range exceeds the redshift interval where the com¬ 
putation of the RS LF is more reliable, i.e. 0.52 < 2 < 0.65. 
Results outside this high-confidence range are obtained by 
extrapolation of the linear fits to the redshift evolution of 
the intrinsic distribution and the fraction of blue objects in 
the sample. 

Figure U displays the redshift evolution of the 3 best-6t 
hyperparameters t, for both the RS and the BC. It is impor¬ 
tant to bear in mind that, due to extreme incompleteness, 
our BC component is not representative of the entire BC 
population, so we simply use it as a means of correcting for 
the contribution of BC objects in the sample. In each panel, 
this BC component is represented in blue lines/symbols. The 
aperture-corrected RS parameters are represented in red and 
the uncorrected RS parameters, in green. 

Statistical errors on each parameter are obtained by 
mapping the likelihood function within the 6-dimension hy¬ 
perspace and computing the posterior probability distribu¬ 
tion for each parameter. This computation yields very small 
errors, which is consistent with the rapid convergence of 
the algorithm. While this procedure might not take all the 
possible sources of error into account, other more realistic 
options, such as a bootstrap analysis, are computationally 
challenging. In order to avoid underestimating our uncer¬ 
tainties, we take a conservative approach here and, instead 
of the aforementioned error, we use, as the real error, the 
scatter on each parameter (standard deviation) with respect 
to the best-fit relation with redshift. 

In the upper panel of Figured the 2.5 -I- ci,_rs increases 
linearly with redshift, from a value of ~ 2.43 at 2 = 0.50 to 
~ 2.53 at 2 = 0.70. The aperture correction has very little 
effect on the zero-point of the logj^Q a-rrii relation. The fol¬ 
lowing linear relation is obtained for the aperture-corrected 
RS Cl (hereafter we drop the ac superindex for simplicity): 

2.5 -k ci,Rs{z) = (2.235 ± 0.009) -k (0.381 ± 0.016) 2 (30) 

At fixed apparent magnitude {mi = 19), we look at 
progressively more luminous galaxies as we move to higher 
redshift. The L-cr relation implies that more luminous RS 
galaxies have higher velocity dispersions, which explains the 
ci.HS - redshift trend. On the other hand, the BC represents 
a photometrically and spectroscopically heterogenous pop¬ 
ulation that contains a large fraction of blue, spiral galaxies 
(or even disky ellipticals) for which we expect smaller veloc¬ 
ity dispersions. At fixed absolute magnitude, the BC com¬ 
ponent is expected, therefore, to have smaller mean logjQ a 
than the RS, which is exactly what the upper panel of Fig¬ 
ure [3] shows. 

The middle panel of Figure [3] shows the redshift evo¬ 
lution of the slope of the logjQ o-mt relation, C 2 . The RS is 
consistent with a single point in the colour-colour plane, with 
a shallow colour-magnitude relation that we neglect in the 


computation of K-corrections. Absolute magnitudes are sim¬ 
ply obtained by rescaling the apparent magnitudes at each 
redshift slice. This result implies that C 2 ,bs coincides with 
the slope of the L-cr relation, as we explicitly show in the 
next section. Our results for C2,rs indicate little evolution 
in the slope of the logj^Q cr-mi within the redshift range con¬ 
sidered. Importantly, the AC, as discussed previously, has a 
significant effect on C 2 ,rs, i.e., that changes from ~ —0.07 
to ~ —0.05. The following linear relation with redshift is 
obtained for the aperture corrected C2,rs'- 


C 2 ,Rs{z) = -(0.033 ± 0.012) - (0.029 ± 0.021) 2 (31) 

As mentioned above, this dependence of C2,rs with red¬ 
shift is not significant given the scatter in the data. The 
Ax^ with respect to the redshift-independent assumption 
is only 1.57 (i.e. a signihcance of 1.25 ct). The slope in the 
logiQ a — mi relation that we measure within the redshift 
range 0.5 < 2 < 0.7 is therefore —0.070 ± 0.006 before AC 
and —0.050 ± 0.007 after. The middle panel of Figure |3] also 
shows in a black dashed line our results for C2,rs assuming 
a constant value for the AC parameter h of b{z) = —0.042. 
Neglecting the redshift evolution of b in Figure |3] has little 
impact on C2,rs- 

The bottom panel of Figure |3] displays the redshift evo¬ 
lution of the intrinsic scatter in the log^Q a — mi relation 
for both the RS and the BC component. This is the in¬ 
trinsic RMS scatter in logj^Q a, at hxed mi. Note that an 
“orthogonal” version of our result (i.e., scatter perpendic¬ 
ular to the logjQ a — mi relation) would be essentially the 
same, given the steep slope that we measure for this relation. 
The scatter for the RS, srs, increases slightly with redshift, 
but this trend is not significant given the computed errors 
(the Ax^ with respect to the redshift-independent assump¬ 
tion is only 1.36, i.e., a significance of 1.17 ct). The mean 
value that we obtain is only srs = 0.047 ± 0.004 in logj^Q cr. 
This value for the scatter in the logj^g a — mi relation is very 
small as compared to the typical value of ~ 0.1 found at 
intermediate-mass ranges and low redshift. In Section |5] we 
show his this result is in excellent agreement with previous 
low-z high-mass results. 

The hyperparameter s for the BC component is signif¬ 
icantly larger, i.e., ~ 0.12. Keep in mind that the BC com¬ 
ponent in our sample is by no means complete, and there¬ 
fore, the corresponding s value that we measure does not 
represent the true scatter for the entire BC population. It 
does provide, however, some indication of the scatter of the 
BC population relative to the RS population. The value of 
~ 0.12 is con s istent with the intrinsic scatter trends shown 
by IShu et ^ (l2012l ~) , who without accounting for the con¬ 
tamination effect caused by the BC objects, find that the 
overall intrinsic scatter value for the CMASS galaxy sample 
is ~ 0.1, an intermediate value due to the mixture of both 
RS and BC galaxies. A trend of h igher intrinsic sc atter at 
higher redshift is also discovered bv IShu et al.l (l2012ll . which 
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Figure 4. The best-fit values for the hyperparameters t that describe the logjQ a-rrii relation for both the RS and the BC component 
as a function of redshift. The distribution of velocity dispersions, for each component, is approximated as a log-normal distribution in 
logj^Q CT with mean < logj^g a > and intrinsic scatter s. The mean is parametrized as < logjg cr >= ci + 2.5 -I- C 2 {rrn — 19). In each panel, 
dots and solid lines represent the best-fit values obtained by maximizing the likelihood function of Equation [21] following our hierarchical 
Bayesian analysis. Blue lines/symbols represent the logj^g a-rrii relation for the BC component, red lines/symbols show this relation tor 
the RS component, once aperture correction is performed, and green lines/symbols show our results before aperture correction. Dashed 
lines show the best-fit linear models obtained within the redshift range considered. The black dashed line in the second panel represents 
a linear fit where aperture correction to C 2 is assumed to be independent of redshift and equal to b(z) = —0.042. Errors on the best-fit 
parameters are assumed to be equal to the standard deviation with respect to the best-fit model in each parameter. 


can be explained by the fact that the relative fraction of the 
BC objects in the CMASS sample increases with redshift. 

The logiQ a-rrii relation in a magnitude-logjg a diagram 
for 4 different redshift slices is explicitly shown in Figure [5] 
(Method 1, see next section). 


7.2 The effect of the various corrections 
implemented 

Our computation of the logj^g a — rrii (equivalent to the L- 
cr relation) at 2 ~ 0.55 incorporates accurate treatments of 
various issues affecting the BOSS data. In particular, we cor¬ 


rect for incompleteness within the colour-colour-magnitude 
space, and we separately model the intrinsic RS and BC 
populations, which allows the statistical removal of all BC 
in the CMASS sample, including those that are predicted to 
scatter through the red side of the colour-colour plane (see 
MD2014 for more details). This information is incorporated 
into our PDSO method, a hierarchical Bayesian statistical 
framework that allows us to utilize the stellar velocity dis¬ 
persion likelihood functions of individual objects, instead of 
the point estimates associated to them. This section is in¬ 
tended to provide a sense as to what impact these corrections 
have on our logj^g a — rrii results. 
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Figure 5. The logj^g a — rrn relation as obtained using 4 different methods in 4 different redshift slices {z = 0.52, 0.55, 0.60, 0.65). Results 
from Method 1, 2, 3, and 4 (see text) are shown in red, blue, black and green, respectively. The solid contours show the best-fit central 
velocity dispersions (point estimates) as a function of apparent magnitude once a colour cut g — i > 2.35 is applied to remove blue 
objects, and the dashed contours represent the entire sample. The best-fit values for the slope and zero-point of the logj^g cr — mi are also 
provided. 


In order to illustrate the effect of each correction, we 
have computed the logjg a — rrii relation using the following 
4 methods: 

• Method 1: This is our optimized method, including all 
the corrections mentioned above. 

• Method 2: This use s exactly the sam e methodology 
as the one implemented in IShu et al.l ll2012l) . applied to the 


dataset used in this work. IShu et ^ (|2012| ) developed of a 
hierarchical approach that incorporates velocity dispersion 
likelihood functions. The main differences with Method 1 
are that 1) the whole CMASS sample is used for the com¬ 
putation, without any red-blue deconvolution, and 2) only a 
very approximative completeness correction is applied. For 
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the sake of comparison, here we exclude blue objects by im¬ 
posing a simple colour cut g — i > 2.35 in observed spaced. 

• Method 3: The observed L-cr relation for the blue sub¬ 
sample dehned using a simple colour cut g — i > 2.35. No 
red-blue deconvolution or completeness correction is applied. 
Instead of velocity dispersion likelihood functions, point es¬ 
timates for the velocity dispersion are used. 

• Method 4: Same as Method 3, but for the entire 
CMASS sample. 


7.2.1 Effect on the slope and the zero-point 

In Figure [S] and Figure |5] we compare in 4 different red- 
shift slices {z = 0.52, 0.55, 0.60, 0.65) the logjQ a-mi relation 
computed using the 4 methods presented above. Figure [5] 
displays these relations in a logj^Q cr vs. rm diagram, where 
the solid contours show the best-fit central velocity disper¬ 
sions (point estimates) as a function of apparent magnitude 
once the colour demarcation is applied, and the dashed con¬ 
tours represent the entire sample. Figure [S] shows the red- 
shift trends for the zero-point and the slope (not aperture- 
corrected, for simplicity), respectively. 

With regard to the zero-point, Figure [5] and especially 
the upper panel of Figure [6l indicate that an inadequate 
or partial removal of blue objects in the sample tends to 
artificially push this parameter towards smaller values. This 
is expected, given that the velocity dispersion is obviously 
smaller in bluer (typically spiral) galaxies (see also the upper 
panel of Figure 3)). A partial removal of blue objects using 
a colour cut (Method 2) only palliates this effect slightly. 
We would still measure a zero point 10% smaller than that 
of the intrinsic RS distribution. Part of this difference can 
also be due to the fact that in Method 2 completeness is also 
only partially addressed, by just applying a rough correction 
to account for the scatter of objects in and out different 
magnitude bins, due to photometric errors. Interestingly, a 
comparison between the zero-point obtained from Method 2 
and Method 3 shows that, for the same sample (both using 
a colour cut in observed space, no deconvolution), the use 
of velocity dispersion likelihood functions in the context of 
a hierarchical Bayesian approach (Method 2) has a minor 
effect on the measured zero-point. 

The effect of corrections on the pre-aperture-corrected 
slope of the logj^g a — mi relation is less obvious. The lower 
panel of Figure |6] shows that the measured slope obtained 
with the optimized Method 1 is only slightly steeper (larger 
in this hgure) than what we would measure using the ob¬ 
served distribution alone (with a colour cut, i.e.. Method 3). 
This differences would likely be within the errors once the 
slope is ape rture corrected ( see following section). Interest¬ 
ingly, using Ishu et al.l ll2012ll method, i.e. Method 2, would 
lead to a pre-aperture-corrected slope much closer to the 
canonical value of 4 (i.e. —0.1 in this hgure), 15 — 25% shal¬ 
lower than the values obtained with Method 1 (note that we 
would still measure a slope steeper than the canonical value 
once aperture correction is applied). 
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Figure 6. Zero-point (upper panel) and slope (lower panel) as 
a function of redshift obtained using the 4 different methods ex¬ 
plained in the text (results from Method 1, 2, 3, and 4 are shown 
in red, blue, black and green, respectively). 


1.2.2 Effect on the scatter 

The use of velocity dispersion likelihood functions in com¬ 
bination with completeness/intrinsic distribution results 
within a hierarchical Bayesian statistical framework has a 
tremendous impact on our ability to recover the intrinsic 
scatter in the logj^g a —mi relation (which coincides with the 
scatter in the L-ct relation). The typical observed scatter in 
Figure [5] ranges from ^ 0.1 at the bright end to ^ 0.16 at the 
faint end. For the ent ire, partia ll y-com pleteness-corrected 
observed distribution, IShu et all (l2012l ) report a value of 
~ 0.1, which is in agreement with the scatter in the higher- 
SN regime (within the sample). Our comprehensive analysis 
allows us to dig even deeper, unveiling the intrinsic scatter 
of the intrinsic RS distribution: a tiny 0.05 dex, as Figure |4] 
shows. 

The difference b etween the scatte r that we measure and 
the one reported bv IShu et al.l ||2012|) illustrates the impor¬ 
tance of the red-blue deconvolution. As shown in the bottom 
panel of Figure|31 the scatter that we expect for the BC pop¬ 
ulation is signihcantly larger. This is consistent with the fact 
that the BC is a much more extended distribution (photo¬ 
metrically and spectroscopically heterogenous) in the colour- 
colour plane. Mixing RS and BC objects and not properly 
correcting for completeness inevitably leads to an increase 
in the reported intrinsic scatter of the L-ct relation. 

In order to demonstrate that we actually have the sta¬ 
tistical power to measure such a small scatter, we have per- 
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formed a simple Monte Carlo variance-estimation analysis to 
estimate the resolution that we can expect given the number 
of objects that we have in a typical redshift bin 20, 000 
objects) and the typical error that we have in onr individnal 
velocity dispersion measurements (~ 0.1 dex, as Figure [T] 
shows). By simulating the observed distribution of veloc¬ 
ity dispersions assuming an intrinsic scatter of 0.047 dex 
in logj^Q (J and the aforementioned typical measurement er¬ 
ror on logjQ a, we can evaluate our capacity to recover the 
intrinsic scatter. Our analysis shows that the typical uncer¬ 
tainty on the intrinsic scatter is of the order of 0.001 dex, 
which is more than one order of magnitude smaller than the 
value that we find for the scatter, i.e. 0.047 dex in logj^g a. 
Note that the error that we estimate for the scatter is 0.004 
dex. This test gives us confidence that our measurement is 
not an artifact that results from working below our resolu¬ 
tion limit. 

7.3 The L-cr relation 

Translating the best-fit hyperparameters shown in Figure [4] 
into the standard L-cr relation (logj^g cr as a function of ab¬ 
solute magnitude) for the RS is straightforward, due to the 
intrinsic characteristics of the RS colour-colour-magnitude 
distribution. In MD2014 we show that, at a given nar¬ 
row redshift slice (width Az = 0.01), and magnitude bin, 
the RS intrinsic distribution is consistent with a single 
point in the colour-colour plane, with only a shallow colour- 
magnitude relation that shifts this point slightly with L. 
Under such conditions, the K-correction, independently of 
the stellar population synthesis model chosen, changes very 
little within the apparent magnitude range of the CMASS 
sample, so we can basically assume it to be constant. We 
can, therefore, convert from apparent magnitudes to abso¬ 
lute magnitudes (K-corrected to z = 0.55) at a given redshift 
slice by simply rescaling the apparent magnitude using the 
standard equation: 

°-^^Mi=mi-DM{z)-°-^^ Ki{z) (32) 

The width of each redshift slice is small enough that 
the variation of DM within the redshift bin can also be 
neglected. By substituting Equation [32] into Equation 1231 
rearranging and adding a factor MqC 2 at each side of the 
equation we arrive at the following expression: 

< logio ^ ^ “f C2 ( — A/q) (^^) 

where: 

c\{z) = ci{z) - C2(z)(19 - DM{z) Ki{z) - Mo) 

C2{z) = C2{z) (34) 

The slope of the L-cr relation is independent of the 
fact that we use apparent magnitudes or absolute magni¬ 
tudes, again due to the characteristic shape of the RS in 


colour-colour-magnitude space. The K-correction as a func¬ 
tion of redshift, ^'^^Ki{z)^ is computed using a grid of mod¬ 
els generated u sing the Flexible St ellar Population Synthe¬ 
sis code IFSPS. fConrov et al.lbOOSlI 'l in the way described in 
MD2014. This grid expands a range of plausible stellar pop¬ 
ulation properties for redshift-dependent models within the 
CMASS redshift range. An average K-correction assuming 
the colours of the RS is computed at every redshift. 

Figure [7] displays parameters cj and C 2 for the L-cr re¬ 
lation as a function of redshift, assuming = —23 

(black). Errors are assumed to be equal to the standard 
deviation of the data for each parameter. The value of 
= —23 for the reference absolute magnitude has been 
chosen because it falls within the CMASS magnitude range 
across the entire redshift range considered (see MD2014). 

As Figure [7] indicates, the zero-point that we measure 
for the L-cr relation has a slight dependence on redshift, so 
that larger values are found at higher redshifts. It is, how¬ 
ever, a very small effect, of ~ 0.005 dex within the red¬ 
shift range considered (0.5 < z < 0.7, ~ 1.3 Gyr of cos¬ 
mic time). The best-fit linear relation that we measure is 
2.429 ± 0.007 -b (0.023 ± 0.011) 2 : (black dashed line). This 
redshift dependence is signihcant as compared to a best- 
fit redshift-independent value of cj = 2.443 ± 0.004 if we 
consider the whole redshift range (Ay;^ = 4.49, i.e., a signif¬ 
icance of ~ 2cr), but not if we restrict ourselves to the high- 
confidence redshift range 0.52 < 2 < 0.65 (Ay;^ = 1.579, 
i.e., a significance of ~ 1.25cr). 

In MD2014 we conclude that the LF evolution of 
the LRG population at 2 ~ 0.55 is consistent with that 
of a passively-evolving population that fades at a rate 
of 1.18 mag per unit redshift. Assuming plausible single- 
stellar population models, in cluding both Flexi ble Stellar 
Population Synthesis (FSPS, IConrov et al.ll2009l 'l and MOO 
(iMaraston et al.ll2009l ~) models, such fading rate, at that red¬ 
shift, translates into a formation redshift for the LRGs of 
2 = 2 — 3. The evolution of the zero-point assuming a best-fit 
passive model of the aforementioned characteristics is shown 
in a red dashed line in Figure (3 Here, by best-fit model we 
mean that the normalization of the zero-point of the pas¬ 
sive model is fit to the data points, so only the redshift 
evolution is meaningful. The deviations found between our 
fiducial model (the best-fit model that we obtained from our 
photometric deconvolution formalism) and the best-fit pas¬ 
sive model are of the order of ±0.004 dex within the entire 
redshift range 0.5 < 2 < 0.7. These deviations are again 
significant (Ax^ = 10.24, i.e., a signihcance of ~ 3.2a) if we 
considered the entire redshift range, but that significance is 
questionable if we restrict ourselves to the high-confidence 
redshift range (Ax^ = 3.01, i.e., a significance of ~ l.73a). 

The fact that the significance of the discrepancies with 
the best-fit passive model reported in MD2014 depend 
strongly on the redshift range considered, in combination 
with the fact the effect is so small (±0.004 dex), appear to 
suggest that our results for the zero-point - redshift trend 
are consistent with the best-fit passive model. This idea is 
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Figure 7. The aperture-corrected L-cr relation parameters as a 
function of redshift within the redshift range 0.5 < z < 0.7. Upper 
panel: The zero-point or, equivalently, the mean of the logj^Q a at 
a K-corrected reference magnitude of = —23 (black dots), 

and the best-fit linear model its redshift evolution (black dashed 
line). In addition, we show the best-fit passive model that fades 
at a rate of 1.18 mag per unit redshift (as measured in MD2014). 
Lower panel: the slope of the L-cr relation, Cj. Errors are assumed 
to be equal to the standard deviation with respect to the best-fit 
model for each parameter. 


reinforced when we consider the uncertainty on the compu¬ 
tation of the AC. A modification of the AC, within reason¬ 
able limits, would produce an effect of the same magnitude 
as the discrepancies that we find between the data and the 
best-fit passive model. As a matter of fact, a variation of 
the a{z) and b{z) functions involved in the AC, within the 
reported uncertainties, could account for the discrepancies 
found. More importantly, the trend shown in Figure[7]is sen¬ 
sitive to the exponent of the AC, as discussed in Section |6l 
We have checked that by increasing slightly this exponent 
we can rapidly reach a much better agreement between the 
model and the data. A value of 0.06 (instead of 0.048), which 
is within the typical range of values previously used in lit¬ 
erature, would produce almost a perfect agreement (a tiny 
Ax^ = 0.11, when using the entire redshift range). 

In summary, we conclude that our results are consistent 
with a passive model with formation redshift z = 2 — 3, 
given the small variation of the zero-point that we measure 
within such a narrow redshift range and the uncertainties 
that the AC is subject to. This result is in good agreement 
with a variety of studies on the evolution of the FP, as we 
discuss in Section EB An improvement in the AC or Re 
estimation in BOSS will be necessary to further constrain 
the evolution of the zero-point. The possibility remains that 


the discrepancies found may also be partially due to the 
fading rate of the passive model being too high; a value closer 
to ~ 0.9 mag per unit redshift would suffice to make these 
discrepancies clearly not significant. This, in any case, would 
not imply a formation redshift much higher than z = 2 — 3, 
given the typical fading-rate evolution of passive models (see 
MD2014 for a discussion). 


8 DISCUSSION 

8.1 Comparison with previous results 

In this section we show that our results for the high-mass L- 
a relation at z ~ 0.55 are in good agreement with previous 
findings at low redshift. 

The canonical form of the L-cr relation for our chosen 
bandpass is Li oc , where Li represents the i-band lumi¬ 
nosity. The slope of this relation, j3, is directly related to the 
measured slope of the a-Mi relation, m^, as 



As we have shown, the mean value of C 2 for the RS, within 
the high-confidence redshift range, is —0.070 ± 0.006, before 
AC is applied. This value implies a L-cr relation slope of 
/3 = 5.71 ± 0.49. Imposing a constant offset of —0.019 to the 
slope, as introduced by the AC according to Equation 1271 
the aperture-corrected L-cr relation becomes even steeper, 
with a mean slope of = 7.83 ± 1.11. 

The above measurement has been performed with a 
sample of more than 600, 000 massive LRGs, with a mean 
stellar mass of M* ~ IO^^’^Mq, within the redshift range 
0.5 < z < 0.7. At intermediate mass ranges and z ~ 0.1, 
using ETC samples extracted from the SDSS, it has been 
shown that the slope of the L-a relation is consistent with 
that of a canonical F-J rela tion, i.e. ~ 4 dBernardi et al.l 
l2003al : iDesroches et al.ll2007l) . Even though this is a fairly 
robust result, studies at these mass ranges might still be sub¬ 
ject to a n inadequate treatment of selection effects. As an 
example. iLa Barbera et~aL I di^) , using the SDSS-UKIDSS 
survey, report a slope of ~ 5 for a sample at a similar mass 
range. 

At higher mass ranges, the curvature of the L-cr relation 
has also been detec ted at high statistical significance at z ~ 
0.1 usi ng the SDSS ([Desroches et al.l[2007l . lHvde fc Bernardil 
l2009al . ISernardi et al.ll201ll ). However, a de finite quantifica¬ 
tion o f the high-mass slope has not emerged. IDesroches et al] 
(|2007^ report a slope of ~ 4.5 at Mr ^ —24 (or logjQ a > 
2.4). The authors find a much steeper slope, of ~ 5.9, for a 
subsample of brightest cluster galaxies (BCGs). 

Some other works have focused on small samples of 
ETCs in the nearby universe (of dozens to a couple of hun¬ 
dred objects). While these samples are strongly affected 
by low-number statistics, they have the advantage that 
galaxy properties can be measured with higher precision. 
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iLauer et alJ ll2007ah . using a compilation of HST observa¬ 
tions of a sample of 219 ETGs, measure a slope for the L-cr 
relation of ~ 7 for a subsample of core and BCG ellipticals, 
the l atter being predomin antly core ellipticals as well (note 
that iBernardi et al.ll200A using an SDSS sample, report a 
significantly steeper slope for BCGs as compare to normal 
ETGs; from their Figure 6 we can visually estimate a slope 
of ~ 8). Core ellipticals are identihed by the fact that the 
central light profile is a shallow power law separated by a 
break from the outer, steep Sersic function profile. A value 
of ~ 2 is reported for the coreless or ’’power-law” ellipti¬ 
cal subsample, a class of objects that present steep cusps in 
surface brightness. Interestingly, a number of studies suggest 
that core ellipticals dominate at the high-mass end, whereas 
coreless ellipticals are predominant at lo wer masses (see e.g . 
iFaber et all 12007 : Lauer et aO 2007lj Hyde et al.l 1200^ . 


Subsequently, iKormendv fc Bender! ll2013ll . using a similar 

sample with some corrections and slight modihcations in the 
core/coreless classification, estimate a slope of ~ 8 for the 
core sample alone, and close to the canonical value of 4 for 
the coreless sample. Our results for the slope of the L-cr rela¬ 
tion are in excellent agreement with these studies, in terms 
of core ellipticals. 


Results consistent with the ab ove picture are also ob - 


tained from the ATLAS®^ sample ( Ca 2 £ellari_et_^ [201l|), 
comprising 260 local-volume ETGs. Cappellari et al. i 2013l l 
measure a high-mass slope of ~ 4.7 (for the related stellar 
mass M-cr relation). This result is obtained above a char¬ 
acteristic stellar mass of ~ 2 x 10^^ Mq, a range of masses 
where the ETG population in this sample is again report¬ 
edly dominated by core ellipticals. Note that a detection of 
this mass s cale in a statistically- siqnificant sample came with 
the SDSS (IBernardi et al.l[201ll L This scale corresponds ap¬ 
proximately to the range of masses covered by BOSS. 


Independently of the central brightness profile clas¬ 
sification, a steep slope at the bright/massive end is 
confir med by other studi e s. In the environment analy¬ 
sis of iFocardi fc Malavasil ll2012l ') (using a small sample 
of a few hundred objects extracted from HYPERLEDA, 
IPaturel et al.l l2003l ~l a value of ~ 5.6 is reported at high 
luminosities, although the authors indicate that an alterna¬ 
tive method to compute this slope could yield a value closer 
to 4.5. 


To summarize, we have drawn two main conclusions in 
terms of the slope of the high-mass L-cr relation from a study 
of previous literature at low redshift. Firstly, enough evi¬ 
dence has been gathered of a steeper slope at the high-mass 
end from SDSS works at z ~ 0.1. Secondly, the reported 
value for the slope varies among different works, from ~ 4.5 
to ~ 8. Note that the L-ct relation is sensitive to low-number 
/ selection effects (which nearby samples are affected by), 
but also to the region of the FP probed (in particular, to the 
exact luminosity range under analysis). Our work, provides, 
for the hrst time, a measurement at an unprecedented statis¬ 
tical significance of the slope at an intermediate redshift and 
at the highest-mass end, where the L-cr relation “saturates” 


(using the terminology of iKormendv fc Benderir2013l . among 
other authors). 

Another important result from this work is the small 
intrinsic scatter of the L-ct relation at the high-mass end 
at 2 ~ 0.55; a measurement that has been performed with 
high statistical signihcance for the first time. The scatter, 
quantihed by the hyperparameter s in this work, is found to 
have a mean value of s = 0.047 ± 0.004 in logj^Q a at hxed L 
and redshift slice (with no signihcant redshift dependence). 
Again, in order to place these results in the context of previ¬ 
ous literature we can only compare with nearby/low-redshift 
samples, where some clear indications have been reported 
that the intrinsic scatter is smaller at higher masses. 

Figure 9 (right panel) from iHvde fc Bernardil (l2009ah 
clearly shows that the scatter in the L-cr relation (in par¬ 
ticular, in the logj^Q ct - Mr relation) decreases towards high 
luminosities. In particular, we can visually estimate a value 
of ~ 0.05 dex at Mr < —23. A similar value for the scatter 
of the mass-CT relation at t he high-mass end can be visually 
estimated from Figure 1 of iBernardi et all 1 I 2 OIIII . Although 
these values are not explicitly reporte d, they are obtained 
from relatively large SDSS samples (the lBernardi et al.ll201ll 
sample contains ~ 18,000 massive ETGs), which indicates 
that these results are statistically significant. 

In nearby samples, however, low number statistics pre¬ 
vent a reliable estimation of the scatter. In any case, indica¬ 
tions have been reported that the intrinsic scatter in smaller 
in core ellipticals, which would be consistent with our 
measu rement. At the high-mass end, IKormendv fc Benderl 
(l2013ll . by adopting measurement errors of 0.1 mag in mag¬ 
nitude and 0.03 in logj^Q ct, report an intrinsic physical scat¬ 
ter of 0.06 in logjQ CT for core galaxies and 0.10 in logj^Q ct for 
core-less galaxies at a given magnitude. 

Even though we must be careful to compare re¬ 
sults at the high-mass end from a quantitative point of 
view, a mass (or luminosity) dependence of the intrin¬ 
sic scatter of the L- ct relation h ave been reported in 
other works, including Isheth et ^ ll2003ll. iDesroches et al. 
l2007L [Nigoche-Netro et al.l ( 2011h and iFocardi fc Malavasi 

( 2 OI 2 I 1 . 

Finally, our result that the evolution of the zero-point 
of the L-ct relation is consistent with a passive model with a 
formation redshift of 2 = 2 — 3 is in good agreement with the 
redshift evolution o f the FP as measured in gal axy clusters 
up to 2 ~ 1 (see e.g . Ivan Dokkum fc Franxlll 996l at 2 = 0.3 9; 

I Kelson et aL 1997 at 2 = 0.58; jyan Dokkum et al.lll99 8l at 


2 = 0.83: Ivan Dokkum et al. 1 1199 a at 2 = 0.83: IWuvts et al 


2004 at 2 = 0.583 and 2 = 0.83; Ivan Dokkum fc Stanford 


2003 : iHolden et ^l2005l at 2 ~ 1.25). More generally, the re¬ 


sult that the high-mass RS population evolves passively from 


a high formation redshift is in agreement with a wide ar- 
rav of analvses (see. e.g.. Wake et all 20061: Gool et al.ll2008l. 

MD2014 for LF results; Maraston et al. 

[2013 for LRG-SED 

evolution results; |Guo et all 20131. 2014 

for LRG-clustering 

results). 

In light of the above comparison, 

it appears that the 
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high-mass end of the L-a relation not only remains un¬ 
changed within the redshift range 0.5 < z < 0.7, but would 
also be consistent with 2 ~ 0 results. 


8.2 Physical interpretation 


The very steep slope of the L-cr relation at the high-mass 
end implies that the interplay between the different pro¬ 
cesses involved in shaping the evolution of RS galaxies at 
different mass ranges is systematically different. As galaxies 
grow, central velocity dispersions do not change as much as 
expected according to the scaling relations at lower masses. 
This result is consistent with the systematic variation of 
the total mass profile as a fun ction of mass in the cen¬ 
tral region of ETGs found by IShu et ^ (l2015h . so that 
more massive galaxies have shallower profiles. One possi¬ 
bility to explain this behavior is that the relative efficiencies 
of gas cooling and feedback in RS galaxies vary at differ¬ 
ent mass scales. Gas cooling permits baryons to condense in 
the central regions of galaxies, and therefore it is believed 
to make the mass distr i bution more centrally concentrated 

fe.g.iGnedin et al]l2004l: Gustafsson et al]l2006l : lAbadi et al.l 


I 2 OI 0 I : Velliscig et al.|[20l4 l. Heating due to dynamical fric¬ 
tion and supernovae (SN)/Active Galactic Nucleus (AGN) 
feedback, in c ontrast, can soften the central density con¬ 
centration fe.g.lNipoti et al.ll2004l:lRomano-Dfaz et al.ll2008l: 


Governato et akF 2010l: Duffy et al. 2010l: Martizzi et al.l 

2 OI 2 I : iDubois et al.ll2013l : IVelliscig et al.l 120141 1. If feedback 


became more efficient in more massive galaxies, that could 
explai n both the results presented here and in IShu et al] 

A more straightforward explanation comes from a sce¬ 
nario where massive and intermediate-mass ETGs have 
a different evolutionary history. From high-resolution im¬ 
ages of a small number of local ETGs, some evidence has 
been gathered that the high-mass end of the RS popula¬ 
tion might be occupied almost exclusively by core ellipti¬ 
cals, whereas corel ess ellipticals domin at e at intermediate 
and lower mass es jLauer et al.l l2007bl B iHvde et al.l l2008l : 
ICappellari et al.ll2013h . This transition occurs at log^Q M, ^ 
11.2, a mass sc ale that was hrs t detec ted with high-statistical 
significance bv \Bernardi et all feOf ill . Even though this clas¬ 
sification arises from the shape of the central surface bright¬ 
ness profile, several studies have shown that this bimodal¬ 
ity, known as the “E-E dichotomy”, extends to a num¬ 
ber of other galaxy properties: core ellipticals have boxy 
isophotes and are slow rotators while coreless ellipticals 


but a few properties, see e.g. 

Kormendv & Bendeil 

1996 

Faber et al.l Il997l: ILauer et al.l 

2007b a: Hvde et al.f 

2008 

Cappellari et al.ll2013l: Kormendv & Bendeill2013l for a com- 


plete discussion). Importantly, the above mass scale has been 
associated, with high statistical si gnificance, with the cur¬ 


vature of the scaling relations (see Hyde fc Bernardill2009al : 
iBernardi et al.ll201lf ). 

This characterization conhrms the ideas of 


iKormendv fc Bender! 1 19961 . who proposed a revision of 
the Hubble Sequence for elliptical galaxies, where isophote 
shape is used as an implicit indicator of velocity anisotropy. 
The general consensus is that the properties of these two 
distinct populations are the consequence of two different 


evolutionary paths. Massive core ellipticals a ppear to be 
formed th rough major dis s ipatio nl ess mergers (iLauer et al.l 

2007^ _ Bernardi_et_^ _ 201 ll : ICappellari et ah 20131 : 

Kormendv fc Bender] l201f ), whereas the less-massive 


coreless ell ipticals seem to have u ndergone a more complex 
evolution (IKormendv et aHl2009l review evidence that they 
are formed in wet mergers with starbursts). 


This paper provides the most precise measurements of 
the L-a relation at the highest mass range ever probed with 
statistical signihcance. Unfortunately, however, BOSS does 
not provide the type of data required to perform a cen¬ 
tral surface brightness profile analysis that can conhrm that 
the sample is dominated by core ellipticals, as previous re¬ 
sults suggest. A complete analysis of this type would an¬ 
swer the question as to whether core ellipticals are fully 
responsible for the steep sl ope. This option is claimed by 
IKormendv fc Bender! (l2013ll by showing that the slope of 
the L-cr relation is significantly steeper for core galaxies even 
in the luminosity region where core and coreless galaxies 
overlap (recall the statistical limitations of the study). In¬ 
dependently of this discussion, the shallow dependence of a 
on galaxy mass (obtained from a L-cr relation very similar 
to wh at we obtained) is reported in IKormendv fc Bender! 
( 2 OI 3 II to be approximately similar to N-body predictions 
( Nipoti et al.l l2003l : iBovlari-Kolchin et al.l I 2 OO 6 I : iHilz et ^ 
I 2 OI 2 II for dissipationless major mergers. 


Previous literature about the “E-E dichotomy” along 
with the si milarities between ou r me asurements and those 
report ed by iLauer et al.l (l2007fJ) and IKormendv fc Bendeil 
(I 2 OI 3 II suggest that the intrinsic high-mass RS distribution 
characterized in MD2014 (the same distribution for which 
we compute the L-cr relation here) can be identified as a 
core-elliptical population. In MD2014, the intrinsic RS dis¬ 
tribution is photometrically deconvolved from photometric 
errors and selection effects. The resulting distribution is so 
narrow in the colour-colour plane that is consistent with 
a delta function, at Hxed magnitude and narrow redshift 
slice (with a shallow colour-magnitude dependence for its 
location). The second component of the intrinsic model, the 
BC, is a more extended distribution, well described by a 
Gaussian function in the colour-colour plane, upon which 
the RS is superimposed. Our BC is dehned as a background 
distribution including everything not belonging to the very- 
pronounced RS (see MD2014 for details). This BC actually 
extends through the red side of the colour-colour plane. This 
characterization could reflect the “E-E dichotomy” on the 
red side of the colour-colour plane, with this intrinsic BC 
being composed by a large fraction of coreless ellipticals, 
for which more scatter on the colour-colour plane is to be 
expected. Follow-up work will be needed to confirm this pic¬ 
ture. 
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The intrinsic scatter of the L-cr relation at z = 0.55 
adds to the result reported in MD2014 that the intrinsic RS 
distribution is extremely concentrated in the colour-colour 
plane as well (at fixed magnitude and narrow redshift slice), 
which is an indication that there is little variability as far as 
the stellar populations are concerned within this population. 
The small scatter in these relations is consistent with the 
idea that this is an aged quiescent population that has seen 
little activity in a long time. 

The above idea is reinforced by the fact that the high- 
mass L-(t relation at 2 = 0.55 appears to be very similar to 
that reported at 2 = 0.1. Also, the evolution of the zero- 
point within the redshift range 0.5 < 2 < 0.7 is consistent 
with that of a passively-evolving population that formed at 
redshift of 2 = 2 —3, in agreement with the LF-evolution re¬ 
sults shown in MD2014 for the same population. The picture 
of the evolution of the high- mass end of the RS is, however, 
not completely clarified yet. iBernardi et al.l (120151') have re¬ 
cently analyzed the high-mass luminosity and stellar mass 
function evolution from the CMASS to the SDSS, reporting 
a puzzling result. The evolution of this functions appears to 
be “impr essively” passive , when K-l-E corrections computed 
from the lMaraston et al.l (l2009l l models are used. However, 
when matched in comoving number- or luminosity-density, 
the SDSS galaxies are less strongly clustered than CMASS 
galaxies, which is obviously inconsistent with a passive evo¬ 
lution scenario. 


9 CONCLUSIONS AND FUTURE 
APPLICATIONS 

We have measured the intrinsic L-a relation for massive, lu¬ 
minous red galaxies within the redshift range 0.5 < 2 < 0.7. 
We achieve unprecedented precision at the high-mass end 
(M* > IO^^Mq) on the measurement of the parameters 
of the L-(t relation by using a sample of 600,000 galaxies 
from the BOSS CMASS sample (SDSS-III). We have decon¬ 
volved the effects of photometric and spectroscopic uncer¬ 
tainties and red-blue galaxy confusion using a novel hierar¬ 
chical Bayesian formalism that is generally applicable to any 
combination of photometric and spectroscopic observables. 
The main conclusions of our analysis can be summarized as 
follows: 

• At 2 ~ 0.55, the passively-evolved L-a relation at M* > 
IO^^Mq appears to be consistent with that at 2 = 0.1. 

• The slope of the 2 = 0.55 L-cr relation at the high-mass 
end is /3 = 7.83 ± 1.11, corresponding to the canonical form 
Li oc . This value confirms, with the highest statistical sig¬ 
nificance ever achieved, the idea of a curved mass-dependent 
L-cr relation. Scaling relations for the most massive LRGs are 
systematically different than the relations defined at lower 
masses. 

• The intrinsic scatter on the L-cr relation is s = 0.047 ± 
0.004 in logj^Q cr at fixed L. This value confirms, with the 


highest statistical significance ever achieved, that the intrin¬ 
sic scatter decreases as a function of mass. 

• We detect no significant evolution in the slope and scat¬ 
ter of the L-ct relation within the redshift range considered. 
Under a single stellar population assumption, the redshift 
evolution of the zero-point is consistent within the errors 
with that of a passively-evolving galaxy population that 
formed at redshift 2 = 2 — 3. This is in agreement with 
the LF-evolution results reported in MD2014 for the same 
population. 

• Our results, in combination with those reported in 
MD2014, provide an accurate description of the high-mass 
end of the red sequence population at 2 ~ 0.55, which is 
characterized in MD2014 as an extremely narrow popula¬ 
tion in the optical colour-colour plane. 

• Our results for the L-cr relation, in the light of previ¬ 
ous literature, suggest that our high-mass RS distribution 
might be identified with the “core-elliptical” galaxy popula¬ 
tion. In light of the ETG dichotomy, the second component 
identified in MD2014, a much more extended distribution 
upon which the RS is superimposed, would contain a signif¬ 
icant fraction of “coreless” ellipticals towards the red side. 
The larger scatter in colour found for this population would 
be consistent with the evolutionary path that has been pro¬ 
posed for coreless ellipticals. 


The above results lead us to consider followup work 
intended to investigate core-coreless elliptical demographics 
in BOSS. This project will require the use of high-resolution 
data. 

The success of our algorithm for the photometric decon¬ 
volution of spectroscopic observables opens a field of future 
applications, as it can be used to constrain the intrinsic dis¬ 
tribution of a variety of spectroscopically-derived quantities 
in BOSS. In the broader picture, the statistical techniques 
developed in this work and in MD2014 lay the foundations 
for galaxy-evolution studies using other current and future 
dark energy surveys, like eBOSS, which are subject to the 
same type of SN limitations and selection effects that we 
face in BOSS. 

The extensive characterization of the high-mass RS pre¬ 
sented in this work and in MD2014 will be used in com¬ 
bination with N-body numerical simulations to investigate 
the intrinsic clustering properties of this galaxy population, 
along with the intrinsic connection between these galaxies 
and the dark matter haloes that they inhabit. The connec¬ 
tion between galaxies and halos will be performed by apply¬ 
ing the techniques of halo occupation distributio ns (HOD: 
e.g.. iBerlind &: Weinber j|2002l : [Zehayi et al.ll200 5lj an d halo 


abundance matching (HA M: e.g.. Vale fc Ostrikerl l2004l : 
iTruiillo-Gomez et al.l 201 ll ~). This is a novel approach as com¬ 
pared with the previous clustering/halo-galaxy-connection 
studies in BOSS, which have focused on the observed galaxy 
distribution and lacked a proper completeness correction. 

Finally, the velocity-dispersion distributions implied by 
the L-cr relation that we have obtained in this work can 
be used in combination with the luminosity-function results 
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of MD2014 to determine the statistical strong gravitational 
lensing cross section of the CMASS sample. This cross sec¬ 
tion can in turn be used to predict and interpret the inci¬ 
dence of spectroscopic ally selected strong l enses within large 
redsh i ft surveys (e.g ., Boltpn__et_aL 20081. Brownstein et al.l 


redsU i tt surveys (e.g ., Ittolton et aLllzUDg, iJaro wnstem et ai.l 
I 2 OI 2 I . lArne.soti et al.il2012l . iBolton et aLri2012l l. and to de- 
rive constraints on cosmological parameter s from the statis¬ 
tics of gravitationally lensed quasars (e.g., lKochand3l 19961 . 
IChae et al.ll200 j . iMitchell et al.lE00Rh . 
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